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Abstract 

Studies on jet substructure have evolved significantly in recent years. Jet 
substructure is essentially determined by QCD radiations and non-perturbative 
effects. Predictions of jet substructure are usually different among Monte Carlo 
event generators, and are governed by the parton shower algorithm implemented. 
For leading logarithmic parton shower, even though one of the core variables is 
the evolution variable, its choice is not unique. We examine evolution variable 
dependence of the jet substructure by developing a parton shower generator that 
interpolates between different evolution variables using a parameter a. Jet shape 
variables and associated jet rates for quark and gluon jets are used to demonstrate 
the a-dependence of the jet substructure. We find angular ordered shower pre¬ 
dicts wider jets, while relative transverse momentum {p±) ordered shower predicts 
narrower jets. This is qualitatively in agreement with the missing phase space of 
p± ordered showers. Such difference can be reduced by tuning other parameters of 
the showering algorithm, especially in the low energy region, while the difference 
tends to increase for high energy jets. 
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1 Introduction 


The determination of different observables related to QCD jets is essential in study¬ 
ing the outcomes of high energy collision experiments. Successful predictions for such 
jet-variables have been achieved by using a combination of perturbative calculations 
at hxed order, parton shower algorithms, matrix-element and parton-shower matching 
algorithms and hadronization models. Study of jet substructure has also evolved sig- 
nihcantly in recent times Bn- Jet substructure techniques are particularly useful in 
identifying the origin of jet(s) in the hard process (5 T6 , and also in removing contam¬ 
ination from pile-up or underlying event [7, 17 22 . 


The discrimination of quark-initiated jets from gluon-initiated ones is an important 
subject involving jet substructure, and has a lot of potential in improving the search 
for new physics. Different methods for quark-gluon tagging have been devised (23}j26], 


with corresponding performance studies 28 30 for the Large Hadron Collider (LHC). 


Theoretical estimates for the performance of such tagging algorithms are primarily car¬ 
ried out with the help of Monte Carlo (MC) simulation tools, such as, Pythia ^ 32 


Herwig 33,34 and Sherpa 35 . Even though qualitative features are in agreement. 


differences in the predictions of the different MC’s have been noted as far as quantita¬ 
tive estimates of the quark-gluon tagger performance is concerned. The primary reason 
for this can be traced back to the fact that the distribution of observables related to 
gluon jets varies signihcantly across the MC’s, while those for the quark jet are largely 
similar. One possible cause of such a feature might be that while tuning the parameters 
of the MC generators, the precise jet data from the Large Electron-Positron Collider 
(LEP) have been crucial, and at leading order in electron-positron collision, the jet data 
is dominantly from quark-initiated processes. As far as the LEP data is concerned. 
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the properly tuned versions of the MC’s have been successful in achieving very good 
agreement with the jet data and are also consistent among each other, even in the 
soft-collinear and the non-perturbative regions. 

Recent studies carried out by the ATLAS and CMS collaborations indicate that the 
data on certain observables related to quark-gluon tagging lies in between the predictions 


of the two MC generators Pythia and Herwig 30,36,1^. Although it might be difficult 


to pinpoint the reason for such differences in the jet substructure observables predicted 
by different generators, understanding the difference between the central components of 
the MC’s can be useful in developing more precise simulation tools. To this end, at a hrst 
order, if we postpone the consideration of the non-perturbative and underlying event 
effects for simplicity, the substructure of a quark or a gluon jet is governed by the pattern 
of QCD radiation, which is controlled by the parton shower algorithm. One of the core 
variables of a parton shower is the evolution variable, different choices for which are 
made in different MC’s. In this study, our aim is to understand the effect of modifying 
the evolution variable and access its impact on jet substructure observables. We also 
ask the question whether certain choice of evolution variables can better reproduce the 
data on quark-gluon tagging observables, as discussed above. 

With this goal in mind, we simulate jet substructure related observables with the 
following generalized evolution variable: 


Ql = |4.-(1 - z)lV, 


( 1 ) 


where, a is treated as a free parameter. For final state radiation, the above variable 
with a = 1 and —1 correspond to the evolution variables employed in PythiaS and 
Herwig++ respectively. In Sec. we provide further details on the framework used to 
implement this evolution variable in our parton shower program. In Sec. we show 
properties of QCD radiations generated by a given and discuss the correlation 
pattern between such radiation properties and the resulting behaviour of one important 

In Sec. 


jet shape observable, Cj 


m 


27 


we show a-dependence of distributions and 


the associated jet rate observable with tuned values of the parton shower parameters. 
We summarize our hndings in Sec. 

2 Formalism 


The evolution variable for the final state radiation of light partons used in our analysis 
is defined in Eq. Q, where z is the momentum fraction of one of the daughter partons 
and is the virtuality of the mother parton. The daughter partons are taken to be 
on-shell. The variable Qa is parametrized by a continuous parameter a, and we take 
the range as a G [—1,1] in this study. Qa with a = 1 and —1 correspond to PythiaS’s 
evolution variable (i.e., relative transverse momentum) and Herwig++’s one respectively. 

When we use Qa as a 


QCD radiations are governed by the DGLAP equation pl9 40 


scale variable, the evolution equation takes on a equivalent form for each a due to the 
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following relation, 




Ql 




( 2 ) 


We implement the general evolution variable Qa for arbitrary a in a parton shower 
program, and calculate jet substructure observables. Even though there are various 
recent parton shower formalisms, e.g., dipole shower in PythiaS 41 or dipole-antenna 


shower in Vincia 42,43 , we use in this study a traditional formalism based on Refs. 34 


44 , which is used in Herwig++. In the following subsection, we describe the modification 


to the formalism in Refs 34,44 required to have a parton shower with arbitrary a. 


2.1 Phase space 

Consider an emission where a mother parton a branches off into light or massless partons 
b and c (a —)■ be). We give an effective mass to the daughter partons to avoid 
singularities in the splitting functions. Then, upper and lower values of the energy 
fraction of one daughter parton z~^ and z'^ are given by 



where is the virtuality of a when b and c are on-shell, and Ea is the energy of a. This 
gives a condition for the allowed region on the energy fraction ze and Qa as 

S •» + - ze), (4) 

^ a ^ max ^ max 

where Qmax and Qmin are the maximal and minimal values for Qa- These are independent 
of a, and given as 


Qmax Qmin 2?7lqg. 


( 5 ) 


Here, z describes not the energy fraction but the light-cone momentum fraction as in 
Refs. 34,44. However, we have explicitly checked that these are approximately the 
same. Hence we use Eq. Q with a substitutions, ze ^ z in the generation of Qa and 
z. The energy of the partons are known at the end of all branchings. So, we set Qmax 
in Eq. Q to the energy of the initial hard scattering process, i.e., ^/s/2 for the first 
branching, and calculate by taking z as the energy fraction for subsequent branchings. 
These choices ensure the required relation = Qi — Qmin — 0) where p_L is the 
spatial component of the relative transverse momentum for each branchings, as defined 
in Ref. 34,44 . 

The allowed phase spaces in the In z — In Qa plane for each choice of a is illustrated 
in Fig. where the parton energy Ea is fixed at 500 GeV. At leading order, the parton 
branchings occur almost uniformly on this plane. The partons start from a high scale 
and evolve to low scale in timelike branchings, and the smaller a is, the larger the phase 
space becomes in the high scale region. So, when the evolution starts from a high scale, 
initial emissions tend to choose a high scale and soft emission for small a. 
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a = +1.0 
a = + 0.5 
or = ±0.0 
a = - 0.5 
or = -1.0 


Figure 1: The allowed phase space in the In z — In Qa plane for each choice of a, with Ea fixed at 500 
GeV. 


2.2 Starting scale 

We consider final states of either a light quark pair (qq) or a gluon pair {gg), with a 
center of mass energy of ^/s, and set the starting scale for the initial partons to their 
energy in the rest frame of the final state, i.e., a/s/2. This is the maximal choice for the 
starting scale, see Eq. ([^. 

Next, consider the sequential branchings a ^ be and b —)■ de, with the scales of the 
branching given by Qa and Qa,b as; 

Ql ~ [4:z{l - 2r)]" X 2z{l - z)El{l - cosOa), (6) 

Ql f, ~ [Azbil - Zb)^ X 2^/1 - Zb)El{l - cos Ob), (7) 


where 6a and 6f, are the angle between b and c, and d and e respectively. The momentum 
fractions for the branchings a ^ be and b ^ de are given by x and Zb, and the energy 
of a and b are Ea and E^ — ZbEa- By imposing the angular ordering 6a > 6h, we get 


Qa,b ^ Q 


4x(l - z) 
Azb{l - Zb) 


-(a+l)/2 
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( 8 ) 

(9) 


The right-hand side in Eq. ([^ can be greater than the previous scale Qa- To avoid this 
wrong of ordering the scale, we set the starting scale of the daughter parton b as 


Ql, = Q„min(l, z[4z(l - z)]-^^+^^/^). (10) 


The angular ordering is ensured by using this starting scale for a = —1. However, 
angular ordered emission is not ensured for a / — 1. Such emissions are vetoed by hand 
as in Pythia6 [M| . 


4 










2.3 Tunable parameters and other modifications 


We use three parameters asi'^z), '^qg^ cind r^ut in our parton shower program. The hrst 
one is the strong coupling constant at the scale of the Z boson mass. We use one loop 
running of as in our code. The argument of as is set to p± = 2“"[z(l — 
thereby including the effects of subleading terms in the splitting functions. The value 
of as is signihcant to the predictions of jet substructure. Larger values of as lead to 
high scale emissions, and jet shape distributions, e.g., the jet mass distribution shift 
to higher value regions. The value of is set to 0.118 in Herwig++, and about 

0.136 — 0.139 for the hnal state radiation in PyhtiaS. The second variable mqg is the 
effective mass of the light partons and gluons to avoid soft-collinear singularities, which 
was introduced in Sec. EH The third one is defined as 


Qcut Qcut 
^cut "Hi ’ 


Qr 


2m. 


qg 


( 11 ) 


where Qcut is a given scale where the evolution terminates. 

We note in passing that, in our analysis, we neglect g ^ qq branchings for simplicity, 
which affect distributions at the NLL order. 


3 Emission property 


Jet shape observables are important in examining the substructure of QCD jets. One 
of th e recen tly studied jet shape observable is the two-point energy correlation function 
, which can be dehned in the rest frame of a parton pair as 




(J) 


27 


45 


o?) = 5: 


EiEj 

i<j &iet 


9ij 

2sin^ 

2 


( 12 ) 


where Ei and Ej are the energies of the particles labeled by i and j in the jet, T^jet is 
the jet energy, and 9ij is the angle between i and j. The sum runs over all distinct pairs 
of particles in the jet. The dominant contribution to this observable comes from the 
hardest emission in the jet, which is also the hrst emission in the jet 46 . Neglecting all 
other emissions except for the hardest one, we get in the soft limit 


In ~ In z J- /? In 


9 

2 sin 2 


(13) 


where and 9 are the smaller energy fraction and the angle of the hardest emission, 
respectively. Evidently from the above equation, studying the properties of the hrst 
emission in the jet on the z — 9 plane will lead to an understanding of the behaviour of 
this jet shape. 

Figs. and 1^ show the emission probability on the \nz — ln(2sin(0/2)) plane for 
quark and gluon jets respectively. The top, center and bottom rows show the results for 
the hrst, second and third emissions. Here, the second and third emissions refer to the 
emissions from the harder of the two partons produced by the hrst and second emissions 
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Figure 2: Emission probability in the Inz — ln(2 sin(0/2)) plane for quark jets. The top, center and 
bottom rows show the results for the first, second and third emissions, respectively. The second and 
third emission refer to the emissions from the harder parton produced by the first and second emissions 
respectively. 


respectively. We find that the equal-probability curves for the hrst emission plots are 
roughly given by the contours described by 


Const. = 


tt -l- 1 


In z -|- In ( 2 sin - ) . 


(14) 


This is because, the evolution variable, in other words, the ordering variable, in Eq. ([^ 
can be written in the soft limit as 


In Qa = —tI— In 2 ; -|- In | 2 sin - ) -|- Const. 


(15) 


It should be mentioned that the small z regions are more favourable due to larger values 
of the strong coupling constant, as- In the case of a = — 1, the evolution variable is 
given by Q-i ~ E x 2sin(6'/2), where E is the energy of the mother parton. So, high 
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Figure 3: Same as Fig. for gluon jets. 


scales also imply larger angles. As mentioned above, the emissions tend to prefer high 
scales and soft emissions for smaller values of a. This is consistent with the results for 
the hrst emission with a = —1 in Figs. [^and|^ 

Clearly, for the jet shape observable in question, we are mostly interested here in 
the first emission in a jet. When we set the jet radius to R = 0.4, such emissions 
are distributed in the region described by ln(2 sin(6'/2)) < —0.9. The hrst emissions 
often fall outside a narrow jet, especially for small a. Also, such emissions tend to 
be vetoed out in the parton shower-matrix element matching algorithms. Therefore, 
it is also important to look into the subsequent emissions. We hnd that the second 
and the third emissions also have a different distribution for each value of a. This 
fact indicates that parton shower algorithms implementing different evolution variables 
would have different predictions for jet substructure. However, the results in Figs. 
and are obtained with the same set of inputs for the tunable parameters described in 
the previous section [[] for all values of a. In the next section, we employ a procedure to 

“ The distributions in Figs. andare obtained with as{mz) = 0.12, ruqg = IGeV, and rcut = 1- 
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a 

as{mz) 

mqg[GeV] 

I’cut 

+1.0 

0.132 

0.94 

1.00 

+0.5 

0.126 

0.90 

1.00 

±0.0 

0.121 

0.84 

1.05 

-0.5 

0.119 

0.83 

1.16 

-1.0 

0.119 

0.85 

1.25 


Table 1: Tuned values of the parton shower parameters for each choice of a, obtained by fitting 
the InCp'^^ and InCp ''^ distributions for quark jets with R = 0.4 with an e+e“ centre of 

mass energy of ^/s = 200GeV. The reference distributions are calculated by using e+e“ —>■ qq events 
generated by Herwig++ 


fit the values of these parameters for each a separately, and show our results with the 
fitted values of the parton shower parameters for completeness. 


4 The a dependence 


4.1 Jet shape distribution 


Jet shape distributions depend on the parameters as{rnz), rugg, and rcut introduced in 
Sec. 2.3 These parameters are determined by performing a £t of the MC predictions 


to experimental data on several jet observables, for which the e+e“ —)■ n jets data from 
LEP are particularly useful. Performing such a fit to the experimental data is, however, 
beyond the scope of the present study as this would require the implementation of a 
hadronization model in our parton shower code. Since the primary goal of this study is to 
examine between difference between parton shower algorithms using different evolution 
variables, as an alternative to real data, we utilize the qq events generated by 

Herwig++ with hadronization switched off as our 

The and distributions have been used to tune the above param¬ 

eters. As mentioned in Sec. the first emission in the jet has a significant effect on 
the jet shape, which can be parametrized by the momentum fraction 2 ; and the angle 9. 
Therefore, two independent distributions contain the necessary information about 
the jet shapes. Here, we use three variables in order to further examine the f3 dependence 
of the QCD jet substructure. 

Throughout this paper, jets are clustered using the generalized kt algorithm for 


collisions using Fast Jet 3.1.1 47 , the distance measure for which is defined as 


^ ^ ^1-cosi?’ 


(16) 


where R is the jet radius parameter, and we use p = —1. 

We hrstly generate events using five choices for the evolution variable, Qi, Qo.s, Q 01 
Q- 0,5 and Q_i at ^/s = 200GeV, where ^/s denotes the center of mass energy in the 

“‘To be specific, we use Herwig++ 2.7.1 with default tune, for the uu and dd parton level final 
states. 













e'^e~ collisions. We calculate and InCj^'^^ distributions with R = 0.4, 

and find the values of the parameters that minimize the variable computed using 
our results and the mock data generated by Herwig++. Theoretical errors are assigned 
using a fiat distribution for each bin. The best fit values of the parameters are shown in 
Table We see that the larger a is, the larger the tuned value of as{mz) becomes. In 
other words, the Pythia8-like case with Qi prefers a higher value of a.s{mz) compared 
to the Herwig-like case with Q-i. This qualitative behaviour is in agreement with the 
actual implementations found in PythiaS and Herwig++. It should be emphasized that 
the outcomes of this tuning procedure do not entirely reflect the Monte Carlo difference 
between PythiaS and Herwig++, as the the parton shower algorithm implemented in 
PythiaS is different from ours. 

In Fig. 1^ the top row shows the fitted results, and hence the distributions are in 
good agreement with Herwig++ predictions. We also obtained the distributions for a fat 
jet (with R = 1.2) and for gluon jets using the fitted values of the parameters shown in 
Table [T] For the same energy, the gluon jet distributions with R = 0.4 are similar for 
each choice of the evolution variable. Small differences appear in the shapes predicted by 
different choices of a for the fat quark and gluon jets {R = 1.2). Fig. shows the same 
distributions as in Fig. with a higher value of the center of mass energy, ^/s = 1000 
GeV. As we can see, the a-dependence of the shapes is found to be higher for higher 
energy jets. 


4.2 Wideness of soft emissions in jets 

The larger the parameter (3 in is, the larger the differences become in Fig. ^ This 
implies that the wideness of the emissions, especially for the hardest emission in the 
jets, is different for each a. This is because, the larger f3 is, the larger the contribution 
to from the emission angle of the hardest emission becomes, which is understood 
from Eq. (13). 

Associated jet rates defined in Ref. directly reveal the wideness of the emissions 
in jets. Associated jets are jets nearby a hard jet, and are defined by two parameters, 
Ra and Ea- Here, Ra is the maximum allowed angle between the momentum directions 
of the hard jet and the associated jet, and Ea is the minimum energy of the associated 
jet^ We set the value to Ea = 20 GeV in this study. 

A high probability for having no associated jet implies that the probability of wide 
emissions occurring around the hard jet is low. Such probabilities have been obtained 
by using PythiaS, PythiaS, and Herwig++, and it has been found that the no associated 


jet probability predicted by Pythia is higher than the one obtained with Herwig++ 38 


The no associated jet probabilities calculated with Qi, Q 0 . 5 , Qo, Q- 0.5 and Q_i are 
shown in Fig. where, the fitted values of the parameters in Table [T] have been used. 
We can see that no associated jet probabilities are similar for each a at the low energy 
range. This is expected as the parameters have been tuned at -^/s = 200 GeV. The 


In Ref. for studies in hadron collisions, the parameter pa has been used to define associated 
jets instead of Ea, where Pa is the minimum transverse momentum of the associated jets. However, for 
the e“'"e“ collisions studied in our paper, it is more suitable to use the energy variable. 
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Figure 4: Distributions of and for quark and gluon jets, with R = 0.4 and 

1.2, at -^/s = 200 GeV, as obtained using the parameter values shown in Table 
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Figure 5: Same as Fig. with a higher center of mass energy, a/s = 1000 GeV. 
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Figure 6: No associated jet probabilities for {R,Ra) = (0.2, 0.4) and (0.4, 0.8), computed with the 
input parameters as in Table 
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Figure 7: Same as Fig. with the input parameters obtained by htting the InC^^'^^ InC^'^^ and 
lnCj^'°^ distributions for quark jets in e+e“ —>■ qq events at i/s = 2000 GeV. 


a dependence is enhanced at the high energy range. The larger a is, the larger the 
no associated jet probabilities become. Therefore, an angular ordered shower (a = — 1) 
predicts wider jets, while apj_ ordered shower (a = 1) predicts narrower jets. This result 


is qualitatively in agreement with the missing phase space of the p± ordered shower 48 


The wideness of the emissions in the jets are thus tunable by changing the parameter a 


in the evo 

Fig. I? 


ution variable continuously. 


is similar to Fig. 


(/3) 


with the tuning parameters obtained by htting InCj 
distributions for quark jets in —)• qq events at y/s = 2000 GeV. The no associated 

jet probabilities are similar for each a around ^/s = 2000 GeV for the quark jets. The 
a dependence now appears at other energy ranges. The energy scaling of the wideness 
seems to be inherent in the choice of the evolution variable for the same modelling of 
the parton shower. 
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5 Conclusions 

In this paper, we have introduced a generalized evolution variable Qa which is a function 
of the free parameter a taking continuous values. Although the evolution equation 
governing the QCD radiation in jets takes an equivalent form for each a, jet substructure 
depends on a even in the same parton shower formalism. We have examined the a- 
dependence of distributions and the associated jet probability for quark and gluon 
jets. This is motivated by the differences found in the prediction for jet substructure 
observables between often-used Monte Carlo generators, and also by the fact that recent 
LHC data related to QCD jet substructure lies between the predictions of the MC 
generators. The angular-ordered parton shower formalism used in this study is built 
upon the one implemented in Herwig++. We leave further studies based on other recent 
parton shower formalisms to a future work. 

We have studied the distributions of the first, second and third emissions in the 
momentum fraction and emission angle 9 plane. These distributions are of importance 
as the beginning emissions in the jets have a signihcant impact on and other jet 
shape observables. The distributions show a unique emission pattern for each choice of 
a. 

We have tuned the parameters in the parton shower to —)■ qq mock data 

generated using Herwig++, with center of mass energies of ^/s = 200 GeV and 2000 GeV. 
Observables used in the tuning are InCj^'^^ lnGp °^ and InCj^'^^ distributions with the 
jet cone angle R = 0.4. From this £t, we observe that larger values of the strong coupling 
are preferred as we vary the values of a from —1 to 1. This is qualitatively in agreement 
with previous findings regarding the difference between the parton shower phase-space 
covered by the p_[_ ordered and angular ordered showering algorithms. Using the best 
fit parameters, we have calculated the InCj^^ distributions of the quark and gluon jets, 
with R = 0.4 and 1.2, for collisions at ^/s = 200 and 1000 GeV. As we move 
away from the setup used for the fits (namely, quark jets, R = 0.4, ^/s = 200 GeV), the 
a-dependence becomes more apparent, especially for larger values of /3 in C^\ 

The a-dependence for large (5 implies that wideness of the soft emissions, especially 
the first ones in a jet are different for each a. We can examine this wideness directly by 
studying the associated jet probability. A high probability for having no associated jet 
simply means that the probability of wide emissions occurring around a hard jet is low. 
We have found that the larger a is, the larger the no associated jet probability becomes. 
This gives us a qualitative understanding of the generator dependence of associated jet 
rates, especially between PythiaS and Herwig++. Our results open up the possibility 
that we might be able to reproduce the wideness of jets observed in real data by varying 
the value of a in the evolution variable continuously. 
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